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Abstract. We aim to find the optimal site locations for a hypothetical network of 1—3 
triangular gravitational- wave telescopes. We define the following A^-telescope figures 
of merit and construct three corresponding metrics: (a) capability of reconstructing 
the signal polarization; (b) accuracy in source localization; and (c) accuracy in 
reconstructing the parameters of a standard binary source. We also define a combined 
metric that takes into account the three figures of merit with practically equal weight. 
After constructing a geomap of possible telescope sites, we give the optimal 2-telescope 
networks for the four figures of merit separately in example cases where the location of 
the first telescope has been predetermined. We found that the optimal site locations 
for a second telescope based on the combined metric form a ±7° annulus at an angular 
distance of ~ 130° from the location of the first telescope. Based on this result 
we conclude that placing the first telescope to Australia provides the most options 
for optimal site selection when extending the network with a second instrument. 
We suggest geographical regions where a potential second and third telescope could 
be placed to get optimal network performance in terms of our figures of merit. 
Additionally, we use a similar approach to find the optimal location and orientation 
for the proposed LIGO-India detector within a five-detector network with Advanced 
LIGO (Hanford), Advanced LIGO (Livingston), Advanced Virgo, and KAGRA. 
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1. Introduction 

Modern gravitational- wave (GW) detectors are state-of-the-art interferometers with arm 
lengths typically on the kilometer scale. The most advanced network of such instruments 
so far consists of the two detectors of the Laser Interferometer Gravitational-wave 
Observatory (LIGO) [I] at Hanford (Washington, USA) and at Livingston (Louisiana, 
USA), both having 4 km long arms, and the Virgo detector at Cascina (Italy) [2j, having 
a 3 km arm length. All three detectors are currently undergoing major upgrades until 
they will come back online around 2014-2017 as a network of second generation GW 
detectors, named as Advanced LIGO [3] and Advanced Virgo [4], respectively. 

Two additional L-shaped detectors are proposed to become operational shortly after 
2018, and will join the network of second generation detectors: the LIGO-India in India 
[5], and the KAGRA detector in Japan [6j, with 4 and 3 km long arms, respectively. 
The KAGRA detector is to be constructed in the Kamioka mine, which already sets 
the location and orientation of this future instrument. However, the site location and 
the orientation of the arms of LIGO-India have not yet been finalized. This provides an 
opportunity to suggest both based on purely scientific figures of merit. 

According to source models, the detection of GWs is expected shortly after 
the second generation GW detectors start taking data with their proposed nominal 
sensitivity [7j. This opens up the era of GW astronomy, where the expected common 
detection of GWs provides an opportunity for observing distant sources and yet 
unexplored physical processes, making GW detectors powerful tools of astronomical 
research, similar to electromagnetic telescopes. To emphasize this major shift in 
the functionality of these kilometer-scale interferometers, we often use the term GW 
telescopes to refer to future GW detectors during and beyond the second generation 
era. Since the results in this paper correspond to such future instruments, we are going 
to use the term GW telescopes to refer to them, and keep the term GW detectors 
for instruments that primarily aim to make the first detection of GWs. In this 
terminology, advanced (second generation) interferometers should only be categorized 
as GW telescopes after they have undergone the expected shift from the first detection 
to the routine detection era. 

Precision GW astronomy requires GW telescopes even beyond the expected 
sensitivities of second generation interferometers, extending also at low frequencies. 
The Einstein Telescope (ET) project [8j is a good example that aims to pioneer novel 
design configurations, construction methods, and operational techniques, to overcome 
the limitations of advanced interferometers in GW astronomy. The final goal of the 
ET project is to provide a baseline for the future construction of a third generation 
GW telescope. This, according to the conceptual document [9], will consist of three 
interferometers (single-interferometer design) or interferometer pairs (xylophone design) 
with 10 km long arms, built underground in a way to form an equilateral triangle shape. 
The ET is aspiring to have a noise level about one order of magnitude below the noise 
floor of Advanced LIGO in its most sensitive range, and extend sensitivity at lower 
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frequencies [10]. This would make the ET capable of observing GW sources at distances 
of even up to several gigaparsecs, depending on the source type [9]. 

GW telescopes are not equally sensitive to GWs coming from different sky 
directions. Beyond the scientific advantages of improved sensitivity and of having 
multiple interferometers at one telescope site, the equilateral triangle geometry makes 
the directional sensitivity of such a GW telescope cylindrically symmetric around the 
symmetry axis of the triangle for both GW polarizations fTT] . In practice, this means 
that all the scientific figures of merit we can define in terms of GW detection will be 
independent from how the triangle is oriented at a given geographical position. For a 
network of such GW telescopes, the corresponding figures of merit will only depend on 
the geographical positions of the individual GW telescopes relative to each other. 

The scientific advantages of the geometrical design proposed by the ET project sets 
a long term model for future GW telescope development in the third-generation era and 
possibly even beyond. Therefore to preserve the generality of our discussions, from now 
on we will use the term Delta telescope (A-telescope) to refer to a future or hypothetical 
GW telescope that consists of three (9(10 km) scale interferometers sharing the same 
site and forming an equilateral triangle shape, similarly to the proposed design for the 
ET. The sketch of a A-telescope is shown in figure [TJ 

Although there are some promising site candidates based on local noise background 
measurements [9], the future geographical position of the first A-telescope is still 
unknown. Beyond the expected local noise background, one can also consider figures 
of merit of a network of GW telescopes when selecting the best site candidates for any 
number of future telescopes. In doing so, we should first consider what type of GW 
telescopes we should take into account in a future network that will include at least one 
A-telescope. 

We consider a first A-telescope that provides an improvement over the second- 
generation GW telescopes in the limiting sensitivity by more than a factor of 10 over 
their entire frequency band [10J. This yields a far better detector performance, and 
different physics to be explored |TT)1 [T2] . As such an instrument becomes operational 
(after 2020 at earliest [9j), the second-generation GW telescopes with less than 10% 
sensitivity and 0.1% explorable volume will become irrelevant as coherent detection 
partners. Therefore we only consider A^-telescope networks here where all the nodes 
are A-telescopes, when exploring the best site locations based on A^-telescope figures of 
merit. 

In case of a single A-telescope, the first figure of merit that we need to take into 
account when choosing the site location from an allowed set of locations is the expected 
noise background caused by local environmental effects. However, when building a 
network of more than one A-telescope, the selection of the first site already severely 
restricts the optimal location for the second (and third, fourth, etc.) telescope, based 
on the A^-telescope figures of merit. By not taking this into account when selecting the 
site for the first telescope, one might get into a situation where the set of site candidates 
for the additional detectors has already narrowed down to options that are sub-optimal 
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0(10 km) 

Figure 1: The sketch of a hypothetical Delta telescope (A-telescope). The instrument 
consists of three identical interferometers (single-interferometer design) or interferometer 
pairs (xylophone design) with (9(10 km) long arms. The opening angle between the arms 
of each interferometer is 60°, and the interferometers are oriented such that they form 
an equilateral triangle shape. The center of the triangle (marked with x) is placed 
to a geographical position represented by [$, A] latitude and longitude coordinates, 
respectively. Note, that this is the baseline design geometry proposed by the ET design 
document [9]. 

in terms of the TV-telescope figures or merit. Therefore even though there has not been 
any published plans yet to construct more than one A-telescopes, it is critical to keep 
the door open for a future extension to TV > 1 A-telescopes and explore the best site 
candidates for a network of multiple A-telescopes. Thus, we emphasize that the results 
based on TV-telescope figures of merit should be taken into account in the site selection 
of the first A-telescope, which we treat as the primary element constructed in a possible 
future network of GW telescopes. 

In this paper, we suggest optimal locations for 1 — 3 A-telescopes from a set of 
possible sites chosen with respect to basic geographic limitations (e.g. availability 
of solid ground; elevation of the site; etc.) and expected noise levels from local 
environmental effects (such as near-coastal microseism and human activity). In the 
analysis, we intentionally disregard geopolitical considerations to allow for best science 
type optimization. The optimization is carried out in terms of three different N- 
telescope figures of merit, and in terms of a combined example metric that takes into 
account the three figures of merit with practically equal weight. In the optimization, 
for simplicity, we assume that all A-telescopes in the network are identical and have the 
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same sensitivity curve. 

Additionally, we use the same approach to find the optimal location and orientation 
for the proposed LIGO-India detector within a five-detector network with Advanced 
LIGO (Hanford), Advanced LIGO (Livingston), Advanced Virgo, and KAGRA. In this 
case, again for simplicity, we assumed that these L-shaped detectors have practically the 
same noise levels. In the frequency band around 100 Hz where these detectors will have 
their lowest level of noise, this approximation is reasonably accurate (e.g. see figure 5. 
in (9j). 

The paper is organized as follows. In section [2] we discuss how the set of acceptable 
site locations were chosen from a geographical map of the world. In section [3] we 
introduce three A^-telescope figures of merit together with the corresponding metrics, 
and a combined metric defined in the three-dimensional metric space. In section HJ we 
give the optimal 2-telescope networks for the four figures of merit separately in example 
cases where the location of the first telescope has been predetermined. We also suggest 
geographical regions in section [4] where the first three telescope in a network should 
be placed to get optimal network performance in terms of our figures of merit, and to 
provide the most options for optimal site selection when extending the network with an 
additional instrument. Section [5] discusses the results of the site location and orientation 
optimization for the future LIGO-India detector. Finally, we summarize and discuss our 
results in section El 

2. Selection of possible telescope sites 

In the process of constructing a set of possible site locations for a future GW telescope, 
we filtered out geographical positions where either the construction of kilometer-scale 
interferometers would be unfeasible, or where we expect to have permanent local noise 
sources that would strongly limit the scientific performance of the telescope during 
observational runs. 

We started our selection process with a cylindrically projected map of the world, 
where the oceanic surfaces have already been filtered out. We removed the polar regions 
from the map due to the harsh weather conditions that would make the construction 
and operation of a GW telescope unfeasible. We also subtracted all the continental lakes 
from our map. Finally, we applied a cut at an elevation limit of 2000 meters above sea 
level, due to the limited accessibility of such places. 

We applied two additional filters that were based on considerations of local noise 
sources. Due to the high level of excess near field seismic noise caused by oceanic waves 
clashing to the continental shelf, we removed a ^ 100 km wide zone along the coastlines 
of continents and islands. To exclude areas that are permanently affected by intense 
human activity that leads to an increased level of gravity gradient background, we used 
a high-resolution image of Earth's city lights from the Visible Earth Catalog of the 
National Aeronautics and Space Administration (NASA) [T3] . By removing the high- 
intensity pixels of the city lights map from our map of site candidates, we excluded both 
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Figure 2: The map of acceptable regions for the construction of a future A-telescope 
(black), together with the excluded geographical areas: oceanic surfaces (white), 
continental lakes (blue), locations at elevations higher than 2000 m above sea level 
(brown), near-coastal regions within ^ 100 kilometers (green), and urban areas (yellow). 
The remaining locations (black) were discretely sampled and were used as a set of 
acceptable telescope sites for our numerical analysis. 

the urban areas and the most heavily used traffic roads. 

The resulting set of acceptable sites, together with the excluded geographical areas 
are visualized in figure El The remnant geographical areas were discretely sampled to 
get a sub-sample of locations where the geographical positions are ~ 300 km away 
from each other. The resulting 1525 locations were treated as the set of acceptable sites 
for A-telescopes. Note, that future analyses with additional constraints and refinements 
might nuance the local selection of acceptable sites, but will not affect the global picture 
we are about to present here. 

3. TV-telescope figures of merit and metrics 

In this section we introduce the figures of merit that we used to evaluate the scientific 
potential of the different configurations of a A-telescope network. Our framework in 
many ways is similar to the ones introduced in [14J and in fT5] . 

As discussed in section [U the figures of merit for a network of A-telescopes are 
independent from the orientation of the individual telescopes. We also pointed out 
that the horizon distances of A-telescopes for the targeted GW source types are in the 
order of many gigaparsecs, in which the distribution of matter can be approximated 
as homogeneous and isotropic. Due to the rotation of Earth and its orbital motion 
around the Sun, the directional sensitivity of a ground-based GW telescope changes 
with time towards the different directions in the galactic frame. By taking into account 
all these effects, as a first approximation, we only need to consider the positions of the 
GW telescopes in the network relative to each other, and use all-sky averaged figures of 
merit to characterize the network. 
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As an alternative, one can configure a network such as to maximize the observational 
time with the highest directional sensitivity towards a certain sky direction. This 
sky direction could be chosen based on known inhomogeneities or anisotropics in the 
mass distribution of the local Universe. However, in this work we optimized the 
network configurations based on all-sky averaged figures of merit, leaving the local mass 
distribution as the scope of a follow-up paper. 

In the following subsections, we are going to introduce three all-sky averaged 
figures of merit together with the corresponding metrics to characterize the different 
configurations of A-telescope networks. Additionally, we are going to discuss why it is 
unnecessary to consider the average directional sensitivity (or reach [H]) of the network 
as a fourth figure of merit. We summarize our figures of merit and metrics at the end 
of this section, in table [TJ 

3.1. Capability of reconstructing the signal polarization 

Interferometric GW telescopes sense the change in the distance between their test masses 
along the lines of the interferometer arms. Therefore in order to model the strain output 
of a telescope for a GW signal arriving from a given sky direction, one must project 
the GW strain to the reference frame set by the interferometer arms. The projection 
can be given as a linear combination of the + and x polarization components of the 
incoming GW signal, where the linear factors, F + and F x (ranging from to 1, and called 
"antenna factors") depend on the sky direction relative to the interferometer arms, and 
on the angle that sets the orientation of the frame in which the two polarizations are 
defined. This frame is called the "polarization frame", while the corresponding angle is 
called the "polarization angle", and conventionally denoted by ^. 

By choosing a ^ polarization angle, one can plot the antenna factors for any ground- 
based interferometers as functions of the different sky directions given by $ latitude and 
A longitude coordinates in the Earth-centered coordinate system. These two-dimensional 
plots are called " antenna patterns" , and they can be constructed for a network of N A- 
telescopes by using the same polarization angle for all interferometers, and calculating 
the root-mean-square average of the + and x antenna factors for the 3A individual 
interferometers: 



Here the F antenna factors can either all correspond to the + or to the x signal 
polarization. Also, for any number of interferometers and for any sky directions 
one can always choose the polarization angle such that F + is maximized, and F x is 
minimized, and therefore F+ > F x . This special polarization frame is called the 
Dominant Polarization Frame (DPF) (see [16] and [T7j for more details). Normally, when 
constructing the DPF for multiple interferometers, one needs to take into account their 
noise levels relative to each other in the frequency band of interest JTT] . However when 
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all interferometers are identical (as in the case of TV A-telescopes) , this consideration is 
not necessary. 

The sensitivity of a network of telescopes towards a chosen direction is usually 
characterized by the root-mean-square average of i^ etwork and j^etwork However, 
depending on the network configuration, the ratio of ^network an( j ^network w i^hin this 
average can have any value from zero to one (note that i^ etwork > ^network - g a } wa y S 
true in the DPF). This ratio towards a given sky direction characterizes the relative 
network sensitivity to the two GW polarizations, and is conventionally called the network 
alignment factor PS] . The value of this factor determines the ratio of the signal-to-noise 
ratios from the two GW polarization components, assuming that in average their sum- 
square energies are the same. 

Having a F£ etwork /^network rat - Q e q Ug j ^o \ j n ^he j^pp me ans that the network is 
equally sensitive to both polarizations of a GW signal incoming from the given direction. 
On the other hand, if the ^network ^network ra ti G is significantly lower than 1 in the 
DPF, that means that the network will be relatively insensitive to the x polarization 
component of the signal (x being defined in the DPF), and thus the network will have 
a limited capability of detecting the x component of a GW signal incoming from the 
given direction. Since in order to reconstruct the polarization of a GW signal, we 
need to be able to detect and reconstruct the energy content of the signal for both 
polarizations, ^network ^network ^ ^ can b e interpreted as optimal in reconstructing the 
signal polarization for the given sky direction, provided that the overall signal strength 
is sufficient. Therefore to characterize a telescope network configuration in terms of its 
capability of reconstructing the signal polarization, we can use the |i^ etwork — ^network 1 2 
quantity averaged over all sky directions: 

1 = (J^ j> \ F T tW ° rk A ) " F* etwork ($, A) f dnj 2 (2) 

where $ and A are the latitude and longitude coordinates of the different sky 
directions defined in the Earth-centered coordinate system, and dQ = cos$rf$rfA is 
the infinitesimal solid angle towards the [$,A] direction. If we want to construct 
an optimal network of telescopes by maximizing its capability of reconstructing the 
signal polarization, we should minimize the all-sky average of the |i^ etwork — ^network 1 2 
difference. As shown by Eq. [2j the I metric is defined in a way that it should be 
maximized instead. 



3.2. Accuracy in source localization 

Transient GW signals that we expect to detect with a future A-telescope will have a 
duration of 0(1 ms — 1000 sec) in the sensitive band. For such signals, the localization 
of the source in the sky can be achieved by triangulation based on the times of detection 
of the incoming signal with multiple interferometers [15] . 

The detection times can only be measured with an accuracy that is directly 
proportional to the signal-to-noise ratio, and thus to one over the strain noise level of the 
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interferometer within the bandwidth of the signal [18] . This uncertainty in the detection 
time measurement can be translated into uncertainty in the source localization for the 
signal. However, when the goal is to compare the accuracy in source localization for 
different configurations of TV identical telescopes, the noise levels of the interferometers 
do not need to be taken into account when constructing a metric. 

Since the triangulation is based on the differences between the detection times, the 
accuracy in source localization also depends on the distances between the telescope sites 
along the line-of-sight to the GW source. If we want to maximize the average accuracy 
of source localization towards all sky directions, we must place the telescopes of the 
network as far away from each other as possible in the three dimensional space. For 
a 2-telescope network, this makes the distance between the two telescopes an obvious 
choice for a metric to be maximized. For a 3-telescope network, it is the area of the 
triangle formed by the three telescopes, that should be maximized (e.g. see Eq.(l) in 
[19]). In order to have the same units for the 2- and 3-telescope metric, we use the 
square of the distance between the telescopes as the metric in the 2-telescope case. In 
both cases, we denote the metric corresponding to the potential of the network in source 
localization as D. 

In case of the LIGO-India tests, where the network consists of TV = 5 telescopes, 
we choose D to be the area of the triangle with the largest area among the triangles 
formed by all possible combinations of three telescopes in the network. Thus, as a 
simplification, we characterize the N = 5 network with the D metric value of its best 
N = 3 sub-network in terms of accuracy in source localization. 

3.3. Accuracy in reconstructing the parameters of a standard binary source 

The most promising candidates of GW emission observable by future GW telescopes are 
binary systems of black holes (BHs) and/or neutron stars (NSs). Based on the models 
of such systems and the proposed sensitivity curve for ET, future A-telescopes should 
be able to detect the GW emission of these binaries up to several or several tens of 
gigaparsecs, depending on the total mass of the system [9]. 

The waveforms of the GW signals from these sources are known with a relatively 
high accuracy from post-Newtonian calculations even in cases when the spins of the 
binary components are not negligible (20, [21]. This makes binary systems of BHs and/or 
NSs standard sources in the sense that they can be used to characterize the capabilities of 
a GW telescope or a network of GW telescopes in terms of distance reach and accuracy 
in source parameter reconstruction. 

In case of a circular binary of spinless compact stars, the typical parameters to be 
reconstructed from the GW waveform are the chirp mass of the system and the distance 
of the source [20] . The alignment of the orbital plane relative to the line-of-sight is 
reconstructed from the relative power in the + and x components of the incoming 
signal [20] , and thus the performance of the telescope network in terms of this is already 
being tested by the / metric (see section 13.11) . 
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If we choose a configuration of TV identical A-telescopes, we can use the Fisher 
matrix method [22] to calculate the chirp mass and source distance reconstruction errors 
for a standard GW source. The source that we chose is an inspiraling binary of spinless 
neutron stars, each having a mass of 1.4M , placed to a distance of 1 Gpc from Earth, 
and aligned such that the orbital plane is always perpendicular to the line-of-sight. For 
this standard source, the sensitivity curves for the single-interferometer (ET-B in [9]) and 
the xylophone designs (ET-C or ET-D in [9]) are practically the same in the frequency 
band of the emitted GW signal, and thus the results are independent from the specific 
sensitivity curve we choose from these three options (we used ET-B in our calculations). 
Also, as we pointed out earlier, a future A-telescope should have no problem detecting 
such a binary source at a 1 Gpc distance, and reconstructing its chirp mass. 

Generally the orbit of an inspiraling binary is elliptical, and the principal axes of 
the ellipse give a preferred polarization basis (20], [23]. However, when the line-of-sight 
is perpendicular to the orbital plane, the components of the base vectors as defined in 
[20] become singular. Since the incoming GW signal is circularly polarized in our case, 
we can choose any two orthogonal vectors as polarization bases perpendicular to the 
line-of-sight, and resolve the singularity. In choosing such a polarization basis, as an 
approximation, we neglected the change in the direction of the line-of-sight throughout 
the time while the signal is within the sensitive band of the A-telescope. 

Using the waveform of the spinless binary in the Newtonian approximation [20], 
we calculated the components of the Fisher information matrix [24] for each single 
interferometer, and summed it up for all the interferometers in the A-telescope sites. 
We found that it is unnecessary to construct a separate metric based on chirp mass 
and source distance reconstruction errors, because the results are the same when we 
compare the different network configurations using the two metrics. Thus we chose a 
metric that is based on chirp mass reconstruction error only. We also found that the 
chirp mass reconstruction errors only slightly depend on the off-diagonal elements of the 
Fisher information matrix (the results with and without the off-diagonal elements are 
within less than 0.5%). Thus for simplicity, we only used the diagonal elements of the 
Fisher information matrix when calculating the chirp mass reconstruction errors. 

The chirp mass reconstruction error, cjm-> depends on the source location relative 
to the positions of the telescope sites ((Jm = ^m^^))- I n order to identify the best 
telescope network configuration, we used the all-sky average of this quantity: 



so that maximization of R should lead to the minimum uncertainty of the 
reconstructed chirp mass averaged to the whole sky. 

3.4- Why we should not care about the directional sensitivity of the network 

One way to characterize the directional sensitivity of a telescope network is to calculate 
the mean square of its combined antenna factors for the + and x polarizations toward 




(3) 
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the different sky directions. This measure will be independent from the polarization 
angle for any geometry of the interferometer arms, and thus independent from the 
orientations of the individual telescopes. By calculating the root-sum-square average of 
this quantity for the whole sky, we get a metric that characterizes the average directional 
sensitivity of the network for the whole sky, a.k.a. the average reach of the network fT4] . 
For a network of TV A-telescopes consisting of 3N interferometers, the metric is defined 
as: 



where F + = F* etwork and F x = F* etwork , both given in the DPF. Due to EqUJ the 
integral in EqJUcan be split to 3N integrals, each corresponding to a single interferometer 
in the network, and the integrals can be calculated individually. The value of X will 
therefore never depend on the configuration of how the TV telescopes in the network are 
arranged, even if the interferometers have a different geometry from the interferometers 
in a A-telescope. Even though by using a different network configuration we get a 
different shape for the combined antenna pattern, the all-sky integral of that surface 
will be the same for all configurations. This practically means, that by choosing a 
different configuration, we might loose sensitivity toward certain directions, but at the 
same time we increase the sensitivity of the network toward other directions. The net 
effect is that the TV-telescope network will be able to cover the same volume of the local 
Universe, which is good enough in our case when we assume that the accessible Universe 
is homogeneous and isotropic. Thus, it is unnecessary to include the X metric in our 
optimization procedure, and we can use only the other three metrics that we defined, 
denoted as /, D ) and R. Nevertheless, for future publications when targeted science is 
considered, a metric characterizing the directional sensitivity of the network, as well as 
the time evolution of cosmic reach will be important. 

3.5. The combined metric 

The three metrics we described so far, /, D, and i?, together define a three dimensional 
metric space. For each TV-telescope configurations, we can calculate the corresponding 
metric values, and identify a spatial point in this space using the metric values as 
coordinates. The goal of the optimization process for the individual metrics is to find 
the configuration of TV telescopes that corresponds to the highest values of /, D ) and 
i?, respectively, from a number of Monte-Carlo trials. This is equivalent to the goal of 
finding the most distant data points in the metric space along the three dimensions, 
separately. From now on we denote these highest values of the individual metrics that 
we actually find in our sample of different TV-telescope configurations as i max , Anax, and 
i?max, respectively. 

The optimal configuration of an TV-telescope network in a combined metric would 
be the one for which all the corresponding metric values are maximal within the sample. 




(4) 
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Table 1: The figures of merit and their corresponding metrics used in the optimization 
of TV-telescope configurations. For detailed descriptions, see section 03 



Figure of merit 


Metric 


Metric description 


Polarization reconstruction 

Source localization (N=2) 
Source localization (N=3) 
Source localization (N=5) 

Source parameter reconstruction 


/ 

D 
D 
D 

R 


_ i 

(h§ |^ etwork ($,A) -^ etwork (^^)| 2 ^) 2 

Square of the distance between the telescopes 
Area of the triangle formed by the telescopes 
Area of the N = 3 sub-network with highest D 


Combined metric 


C 





Unfortunately, due to the competing nature of the three metrics, there is no such 
configuration, thus, we need to construct a combined metric for such a generalized 
optimization process, where the three individual metrics are combined with specific 
weights. 

In our analysis, we did not want to favor any of the three metrics specifically, 
therefore we normalized the metric values with 7 max , Anax, and i? max , respectively. This 
makes all the metric values dimensionless, and somewhat comparable with each other 
in nature. Using this, the combined metric we defined was: 




C=\ 7^ + 7^ + hf- • (5) 



Specific astronomic searches may require different weighting factors, given our initial 
goal, we used equal weights. 



4. Optimal networks of A-telescopes 

We carried out an optimization procedure in three example cases where the location 
of the first A-telescope was predetermined and a second telescope ranges over the set 
of 1524 acceptable sites (see section [2]). In the three example cases the location of the 
first telescope site was chosen to be [$a,Aa] = [48.5°, 18.7°] in Europe (Example A), 
[$b,A b ] = [38.9°, -98.4°] in North America (Example B), and [$ C ,A C ] = [19.6°, 78.0°] 
in India (Example C). 

We then chose a location for the second A-telescope from the set all 1524 acceptable 
sites, and calculated the values of the /, D ) and R metrics (see section Ej) for the 
corresponding 2-telescope network. The resulting set of /, and R values were then 
normalized with the highest value of the three metrics among the sample of 1524 
elements, denoted as J maxj L> max , and i? max , respectively. The J/J max , D/D m8i ^ and 
R/Rm&x ratios (given in percentages) are shown in figure [3] for Example A (Europe), in 
figure [4] for Example B (USA), and in figure [3] for Example C (India). 
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Figure 3: The colormaps of normalized metric values that we get for a network of 
two identical A-telescopes in the case when the first telescope is placed to [$a><\a] = 
[48.5°, 18.7°] in Europe, and the location of the second telescope is chosen from the set 
of acceptable sites described in section El Each circle in the maps corresponds to one 
of the 1524 site locations, while their colors represent the I/I max (a), D/D max (b), or 
R/Rm&x (c) values, given in percentages, where the "max" indices correspond to the 
highest values of the three metrics within the samples of 1524 elements. 

Using the known {/, D, D, i?} max ratios, we calculated the C metric (see 

section l3?5l) for all the 1524 different 2-telescope configurations, and for all three example 
cases. We again normalized the resulting C values with the highest value in the 
sample, C max , and visualized the C/C max values (given in percentages) for the different 
configurations in figure [6] for Example A (Europe), in figure [7] for Example B (USA), 
and in figure [8] for Example C (India). 

Based on the colormaps shown in figure El figure and figure El we found that the 
optimal site locations for a second telescope based on the combined (C) metric always 
form a ±7° annulus at an angular distance of ~ 130° from the location of the first 
telescope. This empirical result provided an opportunity to associate the number of 
acceptable sites within this annulus to all the 1525 acceptable sites that can be used 
as geographical locations for the first element in an optimal 2-telescope network. The 
resulting values, shown in figure [9l measure the number of different optimal 2-telescope 
configurations that can be achieved by using a given location as the base for the first 
A-telescope in the network. As figure [9] shows, placing the first telescope in Australia 
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Figure 4: The colormaps of normalized metric values that we get for a network of 
two identical A-telescopes in the case when the first telescope is placed to [$b,Ab] = 
[38.9°, —98.4°] in North America, and the location of the second telescope is chosen from 
the set of acceptable sites described in section El Each circle in the maps corresponds 
to one of the 1524 site locations, while their colors represent the I/I mQX (a), D/D m8i ^ 
(b), or i?/i? max (c) values, given in percentages, where the "max" indices correspond to 
the highest values of the three metrics within the samples of 1524 elements. 



provides the most options for optimal site selection when extending the network with a 
second instrument. 

The ±7° annulus of optimal sites in the case when the first A-telescope is placed 
to a geographical position of [$, A] = [—35.8°, 146.9°] in Australia is shown in figure fTOl 
We suggest that the site location for the second telescope should be chosen from within 
this annulus if the first A-telescope is placed to [$,A] = [—35.8°, 146.9°]. Note, that 
placing the first A-telescope of the network to any of the acceptable sites within the 
annulus would result with [<£>, A] = [—35.8°, 146.9°] being in the annulus of optimal sites 
for the second A-telescope. 

By setting the location of the first A-telescope in the network to be in Australia, 
and choosing a site for the second A-telescope in the network from the corresponding 
±7° annulus, we can rerun the optimization procedure for a 3-telescope network with 
these two predetermined site locations, similarly to the previous cases of 2-telescope 
networks with the first site locations predetermined. 

As a first example for a 3-telescope network (Example D), we chose one site at 
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Figure 5: The colormaps of normalized metric values that we get for a network of 
two identical A-telescopes in the case when the first telescope is placed to [$c,Ac] = 
[19.6°, 78.0°] in India, and the location of the second telescope is chosen from the set 
of acceptable sites described in section El Each circle in the maps corresponds to one 
of the 1524 site locations, while their colors represent the I/I mQX (a), D/D max (b), or 
R/Rm&x (c) values, given in percentages, where the "max" indices correspond to the 
highest values of the three metrics within the samples of 1524 elements. 



[$di,A D i] = [-35.8°, 146.9°] in Australia, and one site at [$ D 2,A D2 ] = [38.9°, -98.4°] 
in North America. The resulting plots of {/, D ) i?} max ratios, given in 

percentages, are shown in figure dU As shown in the C/C max ratio plot in figure [T2l the 
optimal location for a third telescope in this network would either be in Central Africa 
or in the North-Central region of South America. 

As a second example for a 3-telescope network (Example E), we chose one site at 
[$ei, A E i] = [-26.2°, 139.2°] in central Australia, and one at [$ E2j A E2 ] = [48.5°, 18.7°] in 
Europe. The {/, D ) D ) i?} max ratios, given in percentages, are shown in figure [T3l 

while the C/C max ratios, also given in percentages, are shown in figure [HI The optimal 
locations for the third telescope in this case are in West Argentina and in West-Central 
Africa, based on the C metric. 
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Figure 6: The colormap of C/C max values for 1524 different network configurations 
of two identical A-telescopes, where C max is the highest value of the combined metric 
within the sample of 1524 elements. The location of the first telescope was chosen to 
be at [$, A] = [48.5°, 18.7°] (marked by x), and the location of the second telescope was 
chosen from the set of acceptable sites described in section El 
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Figure 7: The colormap of C/C max values for 1524 different network configurations 
of two identical A-telescopes, where C max is the highest value of the combined metric 
within the sample of 1524 elements. The location of the first telescope was chosen to 
be at [$, A] = [38.9°, —98.4°] (marked by x), and the location of the second telescope 
was chosen from the set of acceptable sites described in section El 
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Figure 8: The colormap of C/C max values for 1524 different network configurations 
of two identical A-telescopes, where C max is the highest value of the combined metric 
within the sample of 1524 elements. The location of the first telescope was chosen to 
be at [$, A] = [19.6°, 78.0°] (marked by x), and the location of the second telescope was 
chosen from the set of acceptable sites described in section [2 
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Figure 9: This colormap shows the number of acceptable sites within the annulus of 
optimal site locations for a second A-telescope in a 2-telescope network, associated to 
the site location of the first telescope in the network. The map is calculated based on 
the 2-telescope combined (C) metric results, where we found that the optimal sites for a 
second telescope form a ±7° annulus at an angular distance of ^ 130° from the location 
of the first telescope. The circles represent the acceptable site locations described in 
section EJ and their color represent the number of optimal site locations associated to 
them. The highest numbers (up to 350) of such optimal sites for a second telescope are 
associated to locations in Australia. 
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Figure 10: This map shows the annulus of the highest number of optimal site locations 
for a second A-telescope in a 2-telescope network, together with the set of all acceptable 
sites described in section [2] ( x symbols) . The annulus is formed by the set of optimal 
sites (represented by small triangles) derived from the combined (C) metric map of a 
2-telescope network where the first telescope is placed to [$,A] = [—35.8°, 146.9°] in 
Australia (represented by a big black triangle). Choosing a location in Australia (and 
the location at [$, A] = [—35.8°, 146.9°] in particular) for the site of the first A-telescope 
provides the most options for optimal site selection when extending the network with a 
second instrument (see figure [9j). In this case, the site location of the second A-telescope 
of the network should be chosen from the elements of the annulus. 

5. The LIGO-India case in the second generation GW detection network 

In this analysis, we used a similar approach as for the A-telescope networks to find 
the optimal location and orientation for the proposed LIGO-India detector within a 
five-detector network with Advanced LIGO (Hanford), Advanced LIGO (Livingston), 
Advanced Virgo, and KAGRA. In this case, we had taken into account in the calculation 
of the / and R metric that the LIGO-India is proposed to be a single L-shaped 
interferometer with 4 km long arms, similarly to the Advanced LIGO interferometers. 
Also, we used the assumption again that the strain sensitivity of all the interferometers 
in the network are practically the same. 

We first created a new set of acceptable sites for LIGO-India by focusing our filtered 
map (see section [2j) to the Indian region and by resampling the map with a higher 
(~ 27 km) resolution. The resulting discrete map of acceptable sites consisted of 5904 
different geographical locations, covering whole India and even extending beyond the 
country. 

We manually chose 3 site locations in India (in the Southern, in the North-Eastern 
and in the North- Western part of the country) , and used each of them as a hypothetical 
site for LIGO-India. The orientation angle of the L-shaped interferometer was defined 



Optimal Networks of Future Gravitational- Wave Telescopes 



19 




-180-150-120 -90 -60 -30 30 
X (deg) 



90 120 150 180 



180-150-120 -90 -60 -30 30 60 90 120 150 180 
I (deg) 



(a) 



(b) 




-180-150-120 -90 -60 -30 30 60 90 120 150 180 
h (deg) 



(c) 

Figure 11: The colormaps of normalized metric values that we get for a network of 
three identical A-telescopes in the case when one telescope is placed to [$di,Adi] = 
[-35.8°, 146.9°] in Australia, and one to [$ D 2,A D2 ] = [38.9°, -98.4°] in North America 
(marked by x symbols). The location of the third telescope is chosen from the set of 
acceptable sites described in section El Each circle in the maps corresponds to one of the 
1523 site locations, while their colors represent the I/I msx (a), D/D mdb ^ (b), or i?/i? max 
(c) values, given in percentages, where the "max" indices correspond to the highest 
values of the three metrics within the samples of 1523 elements. 



as the angle between the East direction and the bisector of the interferometer arms, 
measured counterclockwise. At each site, we calculated the /, R ) D ) and C metrics 
for every orientation angle from 45° to 135° with a resolution of 0.1°. Note, that the 
D metric is always independent from the orientation angle, and that the symmetry 
of projection of an incoming GW to the arms of an L-shaped interferometer makes it 
unnecessary to test any orientation angles below 45° and above 135°. Based on the 
results, the optimal orientation angle for LIGO-India was robustly found to be between 
57.7 — 58.9° depending on the metric that we maximized (+k x 90° due to symmetry 
reasons with k being an integer). 

In the case of site location optimization, instead of randomly choosing the site 
locations for LIGO-India, we carried out our analysis by testing all the 5904 different 
sites. Using each acceptable site in our list as a hypothetical location for LIGO-India, we 
calculated the /, i?, and C metrics for the five-detector network, taking into account 
the known geographical positions and orientations of the Advanced LIGO Hanford and 
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Figure 12: The colormap of C/C max values for 1523 different network configurations 
of three identical A-telescopes, where C max is the highest value of the combined 
metric within the sample of 1523 elements. The location of the first telescope in 
the network were chosen to be at [$di,Adi] = [—35.8°, 146.9°] in Australia, and at 
[$D2, ^02] = [38.9°, —98.4°] in North America (marked by x symbols). The location of 
the third telescope was chosen from the set of acceptable sites described in section El 



Livingston, Advanced Virgo, and KAGRA interferometers. The orientation angle for 
the LIGO-India interferometer was first set to be 45° and colormaps of /// max , D/D max , 
R/Rmax, and C/C max values were produced, where the "max" indices corresponded to 
the highest value of the given metric within the sample of 5904 elements. 

As we found that the /, D, and C metric values are consistently higher in the 
Southern geographical regions, we repeated testing all the 5904 sites using an orientation 
angle of 58.2°, which we had found to be the optimal angle for a LIGO-India placed to 
the Southern part of India, based on the combined metric, C. The results are visualized 
in figure [15] for the three individual metrics, and in figure [16] for the combined metric. 
Note, that testing the different site locations using a 45° orientation angle resulted with 
very similar maps to the ones shown in figure [THl which shows that the optimal site 
locations within India are also very robust to the orientation angle of the interferometer. 

Based on the /, D, and the C metrics, we found that the rule of thumb is that 
the five-detector network of LIGO-India, Advanced LIGO Hanford and Livingston, 
Advanced Virgo, and KAGRA, becomes more and more optimal by placing LIGO- 
India more to the South within the Indian subcontinent. Even though the R metric 
favors the opposite geographical direction for the LIGO-India site, there is only a slight 
dependence of the network performance in terms of this metric. 

As a final result, we chose an example for an optimal configuration for LIGO-India 
to have the site location at [$ = 10.02°, A = 77.76°] and to have an orientation angle of 
58.2° + /cx90° (see figure IT7I). The corresponding i// max , D/D m8iX , R/Rmax, an d C/C max 
ratios (the "max" index corresponding to the highest values of the metrics within the 
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Figure 13: The colormaps of normalized metric values that we get for a network of 
three identical A-telescopes in the case when one telescope is placed to [$ei,Aei] = 
[-26.2°, 139.2°] in central Australia, and one to [$ E 2,A E2 ] = [48.5°, 18.7°] in Europe 
(marked by x symbols). The location of the third telescope is chosen from the set of 
acceptable sites described in section El Each circle in the maps corresponds to one of the 
1523 site locations, while their colors represent the I/I msx (a), D/D mdb ^ (b), or i?/i? max 
(c) values, given in percentages, where the "max" indices correspond to the highest 
values of the three metrics within the samples of 1523 elements. 

sample) as functions of the detector orientation angle are shown in figure [181 Note, that 
the dominant metric in terms of optimizing the orientation angle of the LIGO-India 
detector, is /, i.e. the capability of the network of reconstructing the polarization of 
a detected GW signal. This / metric can suffer even more than a ~ 25% loss if the 
detector is not oriented optimally. 

6. Discussion 

In this section, we summarize the results presented in section [4] and El and suggest 
geographical regions where possible future A-telescopes and the LIGO-India should be 
placed based on the TV-telescope figures of merit we defined in section [31 We also suggest 
an optimal orientation for LIGO-India based on the same figures of merit. 

The figures of merit we used to characterize and compare different configurations of 
TV-telescope networks were the following (see section [3]): (a) capability of reconstructing 
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Figure 14: The colormap of C/C max values for 1523 different network configurations 
of three identical A-telescopes, where C max is the highest value of the combined 
metric within the sample of 1523 elements. The location of the first telescope in the 
network were chosen to be at [^eijAei] = [—26.2°, 139.2°] in central Australia, and at 
[<&E2, ^£2] = [48.5°, 18.7°] in Europe (marked by x symbols). The location of the third 
telescope was chosen from the set of acceptable sites described in section El 

the polarization of a detected GW signal (measured by the I metric); (b) accuracy 
in localizing the GW source in the sky (measured by the D metric); (c) accuracy in 
reconstructing the parameters of a standard binary source (measured by the R metric) . 
We also constructed a fourth metric (denoted by C) that combines the /, D ) and R 
metrics with practically equal weight (see section 13.51) . 

In section [4] we calculated the /, D ) R ) and C maps in three example cases where 
the location of the first A-telescope in a 2-telescope network was set to be in Europe 
(see figure [3] and [6]) , North America (see figure [4] and [7J) , and India (see figure [5] and [8]) , 
respectively. The maps were always normalized with the highest value of the different 
metrics within the sample of 1524 /, R ) and C values. Based on the combined 
metric results, we found that the optimal site locations for the second telescope in the 
network form a ±7° annulus at an angular distance of ^ 130° from the location of the 
first instrument. Using this result we calculated the number of acceptable sites within 
such annuli associated to all the acceptable sites, and found that the highest number 
of optimal sites correspond to locations in Australia (see figure [9]). Thus, we conclude 
that placing the first telescope to Australia provides the most options (see figure [10]) for 
optimal site selection when extending the network with a second instrument. 

The site location for the second L*e/ia-telescope in the network should be chosen 
from the set of acceptable sites within the annulus of optimal sites associated to the 
location of the first site. We chose two example cases of 3-telescope networks where the 
location of one telescope was chosen to be in Australia, and the second telescope was 
placed to North America (Example D; see figure fTTI and fT2l) . and to Europe (Example 
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Figure 15: The colormaps of I/I max (a), D/D mdb ^ (b), and i?/i? max (c) values in case 
of placing the LIGO-India site at different acceptable site locations in the Indian 
subcontinent. The "max" indices correspond to the highest value of the different 
metrics within the samples of 5904 elements. The orientation angle of the L-shaped 
interferometer (being the angle between the East direction and the bisector of the 
interferometer arms, measured counterclockwise) was optimized based on the C metric, 
and was set to be 58.2°. The optimal orientation angle was found to be robust in terms 
of the site location within India. 




E; see figure [T3l and fT4l) . respectively. The optimal sites for a third telescope in such 
network was found to be on Borneo island and in the North-Central region in South 
America for Example D, and in North Canada and in West-Central Africa for Example 
E. 

Using an India map of acceptable telescope sites with ~ 27 km resolution, we tested 
a total of 5904 different site locations for LIGO-India, taking into account the Advanced 
LIGO Hanford and Livingston, Advanced Virgo, and KAGRA as further members of a 
five-detector network (see section [5]). We found the optimal orientation angle of LIGO- 
India to be ~ 58.2° + k x 90° (k being an integer) practically for all site locations 
within India, and based on the combined metric. Here, the orientation angle is defined 
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Figure 16: This figure shows the same results as figure fT5l but for the C metric. 
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Figure 17: A suggested example for an optimal network of second generation GW 
detectors including LIGO-India. In this example the LIGO-India is placed to [$ = 
10.02°, A = 77.76°] with an orientation angle (being the angle between the East 
direction and the bisector of the interferometer arms, measured counterclockwise) of 
58.2° + k x 90°, k being an integer. The geographical positions and orientation angles 
of the other detectors in the network, ([$, A]; a), were set to ([46.4551°, -119.41°]; 171°) 
for Advanced LIGO Hanford, ([30.56°, -90.77°]; 242°) for Advanced LIGO Livingston, 
([43.63°, 10.5°]; 115.6°) for Advanced Virgo, and ([36.42°, 137.3°]; 75°) for KAGRA. 
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Figure 18: The J/J max (blue), D/D max (red), i?/i? max (green), and C/C max (black) ratios 
(the " max" index corresponding to the highest values of the metrics within the sample) 
as functions of the detector orientation angle for a hypothetical LIGO-India detector 
at a geographical location [$ = 10.02°, A = 77.76°]. Note, that the dominant metric 
in terms of optimizing the orientation angle is /, i.e. the capability of the network of 
reconstructing the polarization of a detected GW signal. This I metric can suffer even 
a ~ 25% loss if the detector is not oriented optimally. 

as the angle between the East direction and the bisector of the arms in the L-shaped 
interferometer, measured counterclockwise. 

Also, by using an optimally oriented configuration, we found the Southern-Central 
part of India to be the most optimal region for the site of LIGO-India (probably the 
South-Eastern part of the Karnataka region, or North- Western part of the Tamil Nadu 
region; see figure fT6l) . An example for an optimal location and orientation for LIGO- 
India is shown in Figure dH Comparing the results shown in figure d5] and figure [T8l we 
can also conclude that based on the network figures of merit, setting the orientation angle 
of the LIGO-India detector has greater significance than choosing the site location of the 
detector within India. In the optimization of the site location, the dominant term will 
be the source localization accuracy of the network (D metric), while in the optimization 
of the orientation angle, the dominant term will be the capability of the network of 
reconstructing the polarization of a detected GW signal (/ metric). 

More efforts are required in the future to overcome the limitations of the work 
presented here. Taking a lot more exclusionary criteria into account, possibly with 
continuously variable acceptance instead of the present yes/no type grid, would allow 
constructing a much more detailed and accurate set of acceptable sites. Applying public 
crowdsourcing in setting up a list of exclusionary criteria and/or inaccessible geographic 
areas would result with such a detailed map of possible constructions sites that probably 
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no group of researchers would be able to achieve. 

Additional figures of merits can easily be added to the optimization process 
depending on the aim of various searches from GW transients to the stochastic GW 
background. One such figure of merit could be the expected duty cycle of the network 
or a sub-network of a given number of telescopes, as proposed by [Hj. This figure of 
merit and potentially others could allow not only comparing different configurations 
of the same TV-telescope network (as presented here) but also networks with different 
number of telescopes, leading to optimization of the number of telescopes in the network 
as well. 

The /, D, and R metrics have already been used in similar forms (without all-sky 
averaging) in previous works (see [IB], [IE], and [25]). These metrics can unquestionably 
be used to find the optimal configurations of A-telescope networks for all given TVs, 
however any of them might need more elaboration if one wants to use them to make 
quantitative comparisons between optimal and sub-optimal configurations, or between 
networks of different numbers of telescopes. Also, one can use a different weighing of 
individual metrics when constructing the combined (C) metric depending on the specific 
preferences one would like to follow, which might also be a subject to change with time 
as we move closer to the realization of future GW telescopes. These are all beyond the 
scope of this paper and should be investigated in a future follow-up paper. 

Finally, the optimization of A-telescope networks could be based on regular or 
Markov chain Monte Carlo simulations using the set of acceptable sites and metrics. 
This would allow a more general testing of networks of A > 3 A-telescopes. Also, 
as we pointed out in section El one could take into account the distribution of mass 
in the local Universe (e.g. using galaxy catalogs), and optimize the telescope network 
by maximizing the observational time with the highest directional sensitivity towards 
certain sky directions. We also leave these to be the scope of a future paper. 
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